On the nature of superconductivity in the t-J model of high-T c superconductors 
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Using an extended dynamical cluster approximation we study superconductivity in the two- 
dimensional t-J model. In analogy to the extended dynamical mean field theory, non-local spin 
fluctuations are treated self-consistently with an additional Bosonic bath. To solve the effective 
cluster problem, we employ an extended non-crossing approximation which allows for transitions 
to the symmetry-broken state. At sufficiently low temperatures we find a stable superconducting 
solution with d-wave order parameter. Upon pairing, the exchange energy is lowered, consistent 
with an exchange-based pairing mechanism. 



Introduction. The nature of pairing in high- 
temperature superconductors (HTSC) remains one 
of the most outstanding open problems in materials 
science. After years of active research, we are still far 
from a complete theoretical understanding of the rich 
physics observed in these materials. Early in the history 
of HTSC it was realized that electronic correlations play 
a crucial role not only for the pairing-mechanism in 
HTSC but also for their unusual normal state behavior 
0. Hence, models describing the behavior of itinerant 
correlated electrons, in particular the Hubbard model 
and its strong-coupling limit, the t-J model, were 
proposed to capture the generic physics of HTSC 0, Q . 
However, despite intensive studies, these models remain 
unsolved in two dimensions. 

Conventional (BCS) superconductors are well under- 
stood. Electrons can lower their potential energy by 
forming Cooper-pairs. Chester [4j realized that this re- 
duction in the electronic potential energy actually cor- 
responds to a decrease in the ionic kinetic energy , thus 
providing a clear link between the pairing mechanism and 
the lattice degrees of freedom. In analogy to Chester's 
idea, Scalapino and White suggested an analogous 
relation for the origin of the pairing energy in HTSC 
based on numerical calculations for the t-J model. They 
proposed, if spin correlations are involved in the pair- 
ing mechanism, that the condensation energy in HTSC 
is provided by the exchange energy J uj\ (Si- Sj). Con- 
sistent with this idea, neutron scattering experiments 
in HTSC have shown the evolution of a low-energy 
spin-fluctuation resonance in the superconducting state, 
which can be related to a reduction in exchange energy 
00. 

This analogy to conventional BCS superconductors 
however, has recently been contrasted by optical exper- 
iments in BSCCO y|. While in BCS superconductors 
the kinetic energy slightly increases upon pairing, these 
experiments have shown that pairing in HTSC occurs 
through a reduction of the electronic kinetic energy. In 
order to unify these two seemingly controversial observa- 
tions, it is important to note that in the t-J model, mag- 
netic exchange has mixed electronic kinetic and Coulom- 



bic origins. In addition, even the kinetic term in the t-J 
Hamiltonian has Coulombic contributions due to the re- 
striction of no double-occupancy of a given site. Thus, 
while the Hubbard model with its purely kinetic and 
Coulombic terms can be related to the optical experi- 
ments, the t-J model has a well defined exchange energy 
and hence is the model of choice to study the possibility 
of exchange based pairing. 

Using QMC simulations within the dynamical cluster 
approximation (DCA) , we have shown that pairing in the 
2D Hubbard model indeed happens through a reduction 
of the electronic kinetic energy consistent with the 
experimental observation. To explore the possibility of 
exchange-based pairing, we study in this letter the nature 
of superconductivity in the 2D t-J model. Its Hamilto- 
nian reads in usual notation 



H=-tJ2 



J 



i 



(i) 



where the projected operators c i(T , c itT act in the reduced 
Hilbert-space with empty and doubly occupied sites only 
and rii = Y, a c ia c ia are their corresponding number op- 
erators. Unless otherwise noted, we use the hopping in- 
tegral t for the energy scale and set J = 0.3t. 

The possibility of d x 2_ y 2 pairing in the 2D t-J model 
was indicated in a number of numerical studies of finite 
size systems (for a review see 10] ) . Only recent numerical 
calculations provided evidence for pairing in relati vely 
large systems for physically relevant values of J/t [Ullll2| . 
The question of whether pairing persists in the infinite 
size model however remains open. Thus, simulations for 
the thermodynamic limit are clearly desirable to further 
support the finite size results. 

Simulations within the dynamical mean field theory 
(DMFT) take place in the thermodynamic limit. An 
extended DMFT (EDMFT) study of the pseudogap in 
the t-J model was presented in Ref. 01 . But the lack 
of non-local correlations in the DMFT method inhibits 
transitions to states with non-local order parameters such 
as d-wave superconductivity. The dynamical cluster ap- 
proximation (DCA) extends the DMFT by mapping the 
original lattice model onto a periodic cluster with linear 
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size L c embedded in a self-consistent host. As a result, 
dynamical correlations up to a range ~ L c are treated 
accurately while the longer-ranged physics is described 
on the mean-field level. 

Method. A detailed discussion of the DCA for mod- 
els with local interactions like the Hubbard model was 
given m Ref. [3, Il5| . The first Brillouin zone is divided 
into N c = cells where D is the dimension, with each 
cell represented by its center wave-vector, the cluster K- 
points. The reduction of the AT-site lattice problem to an 
effective A c -site cluster problem is achieved by coarse- 
graining the compact diagrams of the free energy, i.e. 
averaging over the N/N c wave-vectors k within a cell. 
Consequently, the Feynman diagrams for the lattice self- 
energy collapse onto those of an effective cluster model 
embedded in a Fermionic host which accounts for all fluc- 
tuations arising from hopping of electrons between the 
cluster and the rest of the system. 

In the following we outline the modifications necessary 
to treat models with non-local interactions like the t- 
J model. Since the projected operators in the t-J 
Hamiltonian Q do not obey standard Fermionic com- 
mutation relations, the diagrammatic derivation of the 
DCA cannot be applied to the t-J model. However, the 
general idea of the DCA, i.e. the approximation of irre- 
ducible lattice quantities by their corresponding cluster 
quantities can still be applied to the t-J model. An ex- 
tended algorithm for models, such as the t-J model, with 
extended range interactions was formulated in Ref. [lfij 
for the DCA. For the cluster size N c = 1 this alg orithm 
corresponds to the extended DMFA (EDMFA) 16]. Fluc- 
tuations associated with the non-local interactions are 
accounted for by an additional Bosonic bath which re- 
flects the influence of the host on the two-particle level. 
In the case of the t-J model, it represents the fluctuating 
magnetic fields induced by J. 

In the DCA, the effective cluster model is conveniently 
represented in reciprocal cluster space 



K,cr Q 

K,k,<r 

+ X [ J Q q s Q ■ (*Q+<i + *-Q-q 

ker 



-St 



X 



JqSq ■ S-q - -^PqP-q 



h.c. 



(2) 



The operators CjP create (destroy) Fermions on the clus- 
ter with wave- vector K and Sq (pq) are their corre- 
sponding spin- (density-) operators. The sums over k 
and q run over the momenta within a DCA cell about 
the cluster momenta K and Q. In the DCA, non-local 
parameters are coarse-grained, i.e. £k = N c /N ^~J k e K _|_ k 



and Jq = N c /N q ^Q+q are given by the cell-averages 
of the lattice Fourier-transforms ek and J q of t and J. 
The Fermionic and Bosonic hosts are described by the 



auxiliary operators aj^ and 3?q ; with corresponding en- 
ergies Ak and o-> q . When the system becomes supercon- 
ducting, an additional U(l) symmetry breaking term 



(t) 



Hsc = X [ A ka k i a -kT + h - c - 



(3) 



has to be introduced to simulate the existence of paired 
states in the host. 

The solution of the cluster problem given by H c i + 
Hsc yields the cluster Green function G c (K,w), cluster 
spin susceptibility x c (Q, lu) and their corresponding self- 
energies S(K,w) and M(Q,uj). In the DCA, lattice self- 
energies are approximated by their corresponding cluster 
quantities, i.e. the lattice Green function is written as 



G(k, uj) = [ujt - e K+k T 3 - E(K, u)] 



(4) 



To allow for the possibility of superconducting order we 
used the Nambu-Gorkov matrix representation ^3] with 
Pauli-spin matrices Tj. The 11-component of G is given 
by the normal single-electron Green function G(k, z) = 
((c k |.; c^)) z , while the 12-component is the anomalous 
Green function ^(k, z) = ((c k .»; c_^)) z and only finite in 
the superconducting state. The lattice spin susceptibility 
becomes 



x(q,w) 



M(Q,w) + Jq 



(5) 



Self-consistency is imposed by demanding that the clus- 
ter Green function G c and spin-susceptibility Xc are 
equal to the coarse grained, i.e. cell averaged lattice 
Green function 

k 

= [^T -e- K T 3 -S(K,^)-r(K, W )]- 1 (6) 



and lattice spin-susceptibility 



q 

1 



M(Q,u) + J(Q) + 6(Q,w) 



(7) 



respectively, where the sums run over momenta within 
the DCA cell about the cluster momentum K. The DCA 
self-consistent problem is then solved by iteration. In 
this procedure, the auxiliary parameters of the cluster 
model Ak, w q , Ak, V K k and lQ t q are replaced by the 
hybridization functions r(K, uS) and 0(Q, w) which are 
self-consistently determined from the coarse grained lat- 
tice quantities G, Eq. © and x, Eq. J7J. Hence, the 
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particular form of the auxiliary degrees of freedom in 
Eqs. and 10 is irrelevant and thus does not impose 
any further approximation. 

Although the complexity of the problem is highly re- 
duced by the DCA, the solution of the effective quantum 
cluster problem remains non-trivial. The additional cou- 
pling to the Bosonic bath renders many traditional meth- 
ods used to solve the purely Fcrmionic cluster model like 
e.q. QMC difficult to apply. Here we employ an ex- 
tended version of the non-crossing approximation which 
is based on an infinite resummation of certain classes of 
diagrams in a pertubational expansion in the couplings to 
the Fermionic and Bosonic baths. A detailed derivation 
of the NC A equations for the EDMFT quantum impurity 
model (N c — 1) was given in Ref.flil]. The application 
of the NCA to the Fcrmionic DCA cluster problem was 
discussed in Ref.0] and the modifications necessary to 
treat the superconducting state in Ref.jl^l- A detailed 
discussion of the NCA equations for the extended cluster 
problem embedded in the additional Bosonic bath will be 
presented elsewhere. 

Results. We performed simulations for the cluster size 
N c = 4, the smallest possible cluster size that can capture 
pairing with the <i-wave order parameter in the t-J model. 
Thus, this theory should be regarded as the mean-field 
theory for d-wave superconductivity in the t-J model. In 
particular, we choose the conventional set of cluster K- 
points {(0,0),(7r,0),(0,7r),(7r,7r)}. 

To show that our technique is indeed able to describe 
superconductivity within the t-J model, we present in 
Fig. G] our results for the single-particle spectrum at 20% 
doping (S = 0.20) and temperature T = 0.0348i below 
the superconducting phase transition at T c = 0.0377i. 
The density of states N(uj) shown in Fig. develops 
a pseudogap at the chemical potential (ui — 0), when 
the anomalous Green-function F(K,cj) = ((5k|; c_k|)) 
shown in Fig. ^3, becomes finite. The K-dependence 
of the coarse-grained anomalous Green function, i.e. 
F((ir,0),ui) = — F((0, 7r), ui) and zero otherwise (not 
shown), is consistent with a d s 2_ y 2 -symmetry of the or- 
der parameter. Fig. \I]p shows the doping dependence 
of the transition temperature. T C (S) has a maximum 
T rnax _ Q .0377i « 110K (for t = 0.25eV) at the optimal 
doping S « 0.20 and decreases with both increasing and 
decreasing doping. Due to the breakdown of the NCA 
technique at very low temperatures, the phase-diagram 
cannot be extended beyond the region shown and in par- 
ticular no prediction can be made for the critical doping 
beyond which superconductivity vanishes. 

It is interesting to note that without the additional 
coupling to the Bosonic host we were not able to stabilize 
a superconducting solution, at least not above the lowest 
temperatures we can simulate with our technique. This 
shows the importance of a self-consistent treatment of 
the spin- fluctuations between the cluster and host. 

To explore the nature of pairing in our simula- 




FIG. 1: (A) Density of states and (B) real part of the coarse- 
grained anomalous Green function in the superconducting 
state of the 2D t-J model at 20% doping at T = 0.0348t. 
(C) Doping dependence of the critical temperature T c . 



tions, we calculate the exchange energy E xc from the 
spin-susceptibility x(q, w), Eq. J^Jl, via the fluctuation- 
dissipation theorem as 



E~ 



l£j( q )<£( q ).£(-q)) 
q 

-ilVj(q) fdJ*™*^ 



(8) 



To obtain results for the normal state below the tran- 
sition temperature T c we suppress superconductivity by 
not allowing for a finite order parameter in our simula- 
tion. 

The exchange energies in the normal and superconduct- 
ing states as a function of temperature close to the critical 
temperature T c at 20% doping are shown in Fig. [21 Our 
results for the nature of pairing in the 2D t-J model are 
consistent with an exchange-based pairing mechanism: 
Below T c , the exchange energy in the superconducting 
state is reduced compared to the normal state, while 
the energy coming from the kinetic (first) term in the 
t-J model, Eq. slightly increases upon pairing (not 
shown). It is important to remember that in the map- 
ping of the Hubbard model to the t-J model, the exchange 
interaction J comes from virtual hopping processes be- 
tween neighboring sites and thus has kinetic character. 
Hence the present results are consistent with our results 
for the Hubbard model where we found the supercon- 
ducting instability to be accompanied by a lowering of 
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FIG. 2: Exchange energy as a function of temperature at 20% 
doping across the superconducting transition in the normal 
(NS) and superconducting (SC) states. 



the electronic kinetic energy [9| . The nature of pairing 
in the t-J model however is more similar to that in BCS 
superconductors: The energy of the exchange coupling 
which plays a similar role to that of the lattice degrees 
of freedom in BCS superconductors, is reduced. 

Neutron scattering experiments relate the gain in ex- 
change energy to the development of a resonance in the 
spin-susceptibility near Q = (it, w) and 40 meV. Accord- 
ing to Eq. (SJ, the exchange energy E xc is reduced when 
the spectral weight in the spin-susceptibility x(q, to) is 
increased near q = (tt,tt) since J(q) = J(cosq x + 
cosq y ) < in this region in q-space. In Fig. |31we plot 
the imaginary part of the dynamic coarse-grained spin- 
susceptibility at Q = (n, 7r) in the normal and supercon- 
ducting states at 20% doping for two different tempera- 
tures below the critical temperature T c . This q-space av- 
eraged quantity allows us to infer the changes in spectral 
weight which according to Eq. JSJl reduce the exchange 
energy. Although our data does not display the develop- 
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FIG. 3: Imaginary part of the coarse-grained dynamic spin- 
susceptibility in the 2D t-J model for Q = (n, n) at T=0.036t 
(top) and T=0.034t (bottom) in the superconducting (SC) 
and normal (NS) states. 

ment of the narrow resonance peak at ~ 40 meV observed 



in experiments, we find that the reduction in exchange 
energy indeed originates in an increase of spectral weight 
in Xs(q? w ) m the region near Q = (n, n) consistent with 
experiments. 

To summarize, we have presented extended DCA/NCA 
simulations for the 2D t-J model which show evidence for 
a ci-wave superconducting state in the thermodynamic 
limit. We find that the pairing instability is accompa- 
nied by a reduction of the exchange energy consistent 
with exchange-based pairing. In agreement with neutron 
scattering experiments, this exchange energy reduction 
originates in an increase in the magnetic response around 
the antiferromagnetic wave- vector Q = (tt, it). 
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